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The luminosity function of Lya emitters at2.3<z<4.6 from 
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ABSTRACT 

We have used Visible MultiObject Spectrograph Integral-Field Unit (VIMOS-IFU) obser- 
vations centred on a radio galaxy at z = 2.9 to search for Lya emitters within a comov- 
ing volume of ~ 10 4 Mpc 3 . We find 14 Lya emitters with flux > 1.4 x 10~ 20 W m~ 2 , 
yielding a comoving space density of 0.0018^q'qqq5 Mpc -3 . We fit a Schechter luminos- 
ity function which agrees well with previous studies both at similar redshift (z ~ 3.4) and 
higher redshift (z ~ 5.7). We therefore find no evidence for evolution in the properties of 
Lya emitters between 3 < z < 6, although our sample is small. By summing the star- 
formation rates of the individual Lya emitters we find a total cosmic star-formation rate den- 
sity of psfr = 6.7 ± 0.5 x 10 -3 M0yr _1 Mpc -3 . Integrating over the luminosity function 
for the combined Lya surveys at z ~ 3.4 and accounting for the difference in obscuration 
between the Lya line and the UV-continuum yields an estimate of psfr. ~ 2.2 x 10~ 2 
MQyr -1 Mpc -3 , in line with previous multi-colour and narrow-band surveys of high-redshift 
star-forming galaxies. 

The detection of high-redshift emission-line galaxies in our volumetric search shows that 
the unique capabilities of wide-field integral-field spectroscopy are well suited in searching 
for high-redshift galaxies in a relatively unbiased manner. 

Key words: Cosmology:observations - galaxies:Distances and Redshifts - galaxies:evolution 
- galaxies:formation 



1 INTRODUCTION 



One of the most important questions in cosmology is how the Uni- 
verse evolved into what we see today. It is believed that structure 
forms from the initial density perturbations in the dark matter distri- 
bution which collapse under gravitatio n al instability lEggen et al , 
196 1 ISandage et all 1 197(1 Peebles! 1 197 it IPress & Schechtei . 
1974 JWhite & Reesl 19781) . The gravitational attraction of the dark 



matter overdensity ca uses gas to accumula te and cool, thus en- 
abling star formation lWhite&Reeslll978l) . Hierarchical models 
of galaxy formation postulate that these sub-galactic clumps subse- 
quently merge to form larger galaxies. Studying star-forming galax- 
ies as a function of both luminosity (or star-formation rate) and 
redshift is therefore essential to obtain insight into the formation 
history of galaxies. 

As galaxies form and evolve they also influence their envi- 
ronment. The most obvious way in which this happens is that 
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they are able to ionise the gas surrounding them. Observations 
of distant quasi stellar objects (QSOs) exhibit absorption short- 
ward of the Lya line iGunn & Petersonf Il965t) . which indicates 
that the reionisation of the Universe was not yet c omplete at z > 6 
(e.g. lBecker et all 1200 lllDiorgovski et all 1200 ll) . The vital ques- 
tion remains: which objects are responsible for this process? Again 
analysing the properties of primeval galaxies can provide an answer 
as the radiation emitted by the young OB-stars coul d be energetic 
enough to ionise the intergalactic medium (see e.g. lMadau et all 

To detect primeval galaxies and determine the cosmic star- 
formation rate density, two different types of object are often tar- 
geted. One comprises the Lyman Break Galaxies, which are star- 
forming galaxies characterised by a strong break in their spectrum 
shortward of the Lyman limit at 912A and the hydrogen Lya line 
at 1216A. The other type is the focus of this paper: the Lya emit- 
ters, identified by their prominent emission in the Lya line. Lyman 
Break Galaxies have been discovered in hu ge numbers over th e 
past decade using multi-colour surveys (e.g. lSteidel etall Il999h . 
Using this technique several surveys, probing different ranges in 
redshift, indicate a rou ghly constant star-formation rate density be- 
tween 1 < z < 5 Iconnollv etail 1 19971 : iMadau et all 1 19981 
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ISteidel et all Il999t llwata et all l2003h . Ho wever, at z > 6 the 
star-formatio n rate density seems to decline jStanwav et afll2003t 
iBunker et alll2004 possibly indicating that the source of reionisa- 
tion may not be the star-forming galaxies we see at 2 < 6. Searches 
for Lya emitters are complementary in this goal as these surveys 
can probe galaxies with relatively very weak continuum emission - 
undetected by the Lyman Break technique. 

Several techniques have been used to observe high-redshif t 
Lya emitters: (i) slitless spectroscopy (e .g. [Koo^^Crorj. ^8C|). 
(ii) sl it-spectroscopy (e.g. lLowenthal et"anft^97I ^hom£son^taU 
1 19921). and (iii) n arrow-band imaging (e.g. iHu et al.(. 1 19971 : 
ICowie&Hulll998t) . Each method has its pros and cons. Slitless 
spectroscopy allows a large volume to be probed, but needs very 
long integration times to achieve an acceptable flux limit. Slit- 
spectroscopy can achieve a good flux limit with shorter integration 
times but unfortunately the surveyed volume is often very small 
and multi-colour imaging observations are usually needed to pre- 
select candidates. Narrow-band imaging compares broad-band im- 
ages to images taken with a filter of very narrow bandwidth (e.g. 
AA ~ 60A). Objects that are significantly brighter in the narrow- 
band image are likely to have a strong emission line at the observed 
wavelength of the narrow-band filter. This way large areas can be 
searched for Lya emitters. This type of survey has been the most 
successful in detecting Lya emitters since the late 1990s when the 
sensitivity of the surveys increased considerably due to the on- 
set of 8-metre and 10-metre class telescopes (e.g. IHu et alj.ll997t 
ICowie & Hulll998h . 

With this leap in collecting area, detection of Lya emitters 
through narrow-band imaging has frequently been used to deter- 
mine the cosmi c star-formation rate densit y from z m 3 — 7 (e.g . 
iKudritzkietail l200fi iFuiita et all l2003t iKodaira et all 120031) . 
Nevertheless narrow-band imaging suffers from one major draw- 
back: due to the small FWHM of the narrow-band filter, the surveys 
sample a very narrow redshift range, and thus a shallow volume. 
Hence it is impossible to use a single narrow-band survey to deter- 
mine the Lya density and luminosity function throughout a large 
redshift range. 

In this paper we aim to improve upon these limitations by 
searching for Lya emitters with integral-field spectroscopy. An 
Integral-Field Unit (IFU) samples an area on the sky where every 
pixel contains a spectrum. This then allows a volume of the Uni- 
verse to be probed with just a single observation. Thus depending 
on the wavelength range probed we are able to study the comoving 
space density and properties of Lya emitters over a large redshift 
range. 

We have used the Vi sible MultiObject Spe ctrograph Integral- 
Field Unit (VIMOS-IFU, LeFevre et aUl2003h on the Very Large 
Telescope to search for Lya emitters within a volume of ~ 10 4 
Mpc 3 . In section 2 we describe our observations and the data re- 
duction and section 3 details the selection procedure used to find 
Lya emitters in our data. The results are presented in section 4; 
section 5 discusses our findings and in section 6 we provide a 
summary of our conclusions. Throughout this paper we assume 
^matter = 0.3, J7a = 0.7, and Ho — 71 km s _1 Mpc _1 , unless 
stated otherwise. 



2 OBSERVATIONS AND DATA REDUCTION 

All of our observations were made on the nights of 29 April 2003 to 
2 May 2003 on the 8.2-metre Melipal telescope at the ESO Paranal 
Observatory in Chile. We used the VIMOS-IFU which consists of 



13" 




j;: 1 « 


!=■ | ■ | | m | | ■ | | 


28" 






dur t i!K , 


C 43" 
O 






4 1-0 










O 








c 








(U 28' 58" 
Q 








13" 














4 , ■ . 


-24° 29' 28" 









09 h 45 m 35.4 s 34.4 s 33.4 s 32.4 s 31.4 s 30.5 s 

Right Ascension 
Figure 1. The two-dimensional image of our IFU field of view, integrated 
from 4050 A to 6800 A. 



6400 fibres coupled to microlenses. Each fibre produces a spec- 
trum of a small region on the sky; these spectra can be integrated 
to construct a 'broad-band' image, as shown in Fig.Q Our obser- 
vations were made in low-resolution mode which corresponds to a 
resolution of R ~ 250. The spectral wavelength range was 3500A 
— > 7000A with a dispersion of 5.35 A per pixel. The spatial sam- 
pling of the fibres in our set-up was 0.67 arcsec per fibre, resulting 
in a square continuous field of view of 54 x 54 arcsec 2 with the 
dead space between the fibres comprising less than 10% of the fi- 
bre to fibre distance. The field is split up into four quadrants, each 
connected with an EEV CCD detector consisting of 2048 x 4096 
pixels. 

The observations were centred on MRC 0943-242, a powerful 
radio galaxy at z = 2.92 (a = 9 h 45 m 32.7 s , S = -24°28'49.7", 
J2000). The total exposure time was 9 hours, split into 18 expo- 
sures of 30 minutes. In order to minimise the effects of bad pixels 
and contamination by cosmic rays during data reduction, the con- 
secutive exposures were offset by 10 arcsec in both right ascension 
and declination relative to the centre. During all four nights the see- 
ing varied between 0.6 and 0.8 arcsec. 

The data were reduced using the V IMOS Interactive Pipeline 
Graphical Interface (VIPGI, IScodeggio et alll2004h . First the stan- 
dard routines for bias and flat-field correction are applied. Subse- 
quently the individual spectra on the science frames are separated 
from each other. Each frame consists of 400 spectra, correspond- 
ing to one quadrant of the IFU. The layout of the IFU fibres in the 
instrument focal plane is used to calculate the position of the in- 
dividual spectra on the CCD frame. When the spectra are located, 
the wavelength calibration is derived specifically for each fibre by 
identifying arc lamp lines taken with the same set-up. The varia- 
tions in the transmission of the individual fibres are corrected for 
and a sky estimate is calculated by averaging the individual spec- 
tra. The sky spectrum is then subtracted from the individual spectra 
and cosmic rays are removed using a 3a clipping algorithm on the 
data. The spectra are flux calibrated using standard-star frames. To 
this end the instrument sensitivity function is derived by summing 
the spectra which contain flux from the standard star. This results 
in an accurate relative flux calibration. To obtain the absolute flux 
calibration, we use a long-slit spectrum of the central radio galaxy 
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taken with the LRIS on the Keck Telescope. We measure the flux 
of the Lya line by integrating the flux in the LRIS-spectrum and 
construct a spectrum of the same object with the VIMOS-IFU data 
by summing the fibres along the long-slit orientation. This allows 
us to bootstrap the spectrophotometry. Comparisons of spectra and 
V-band photometry of several objects in the field show that the fi- 
nal flux calibration is good to ~ 20 per cent. The effect of the dead 
space between the fibres on the error in the flux measurement is 
negligible as it is inherently small and has been even further re- 
duced by the dithering of our single exposures. 

The science frames are median combined to remove any re- 
maining cosmic rays and a fringe correction is applied by median 
combining the fringe frames. The final data cube consists of 12100 
combined spectra, as a result from our enlarged field of view due 
to the dithering of the individual exposures. These 12100 spectra 
correspond to a two-dimensional image of 110 x 110 pixels with a 
field of view of 1.2' x 1.2' (see Fig.0. 



3 DATA ANALYSIS 

3.1 The sensitivity of the data 

To determine the sensitivity of the data a wavelength range is se- 
lected where the flux calibration is good and the spectra are clear 
of skylines. At both ends of the spectra the flux calibration is poor, 
so we discard the regions at 3500A < A < 4200A and 6800A 
< A < 7000A. Between A = 5500A and A = 6550A three sky- 
lines dominate the spectra; this region is also omitted from the de- 
termination of the sensitivity function. In each spectrum we fit a 
polynomial to the usable wavelength range to facilitate calculation 
of the Root Mean Square (RMS) of the noise (a). The fit is sub- 
tracted from the data and a is calculated. As the presence of actual 
features in the spectrum would alter the value of the noise, data 
points deviating more then 3a from the fit are discarded. Subse- 
quently the RMS of the noise is recalculated. This procedure results 
in a measurement of the overall RMS of the noise in each spectrum. 

Visual inspection of random spectra with the continuum fit 
and noise estimate confirm that the measurements are accurate. 
The method does not work well in fibres where there are bright ob- 
jects, e.g. stars or nearby galaxies, due to the large number of emis- 
sion and/or absorption lines in the spectrum. However, to search 
for distant Lya emitters we want to avoid those regions of the 
sky with bright nearby objects. Therefore the sensitivity determi- 
nation described above is well suited to our needs as it yields an 
accurate value for the noise in spectra that do not sample lumi- 
nous sources. To be able to exclude regions with bright objects the 
noise value in the corresponding fibres is set to zero. Inspection 
of the spectra of these objects reveals them to be nearby galax- 
ies, st ars, a QSO at z = 1.8 (see ljarvis. van Breukelen. & Wilmarl 
120051) and a galaxy at z = 0.62, identified through its [Oil] emis- 
sion. In Fig. [2] the sensitivity function over the two-dimensional 
image is shown. Note that, especially in the smoothed image, the 
effects of dithering the exposures are evident. The noise varies from 
a = 1.4 x 10""" 
the data cube. 
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3.2 Selection of line emitters 

With a dispersion of 5. 35 A per pixel and a resolution of R m 250, 
the FWHM of the emission lines will generally be only about 3 
pixels as we expect the typical Lya emitter to have an unresolved 
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Figure 2. The sensitivity function over our field of view. The lighter re- 
gions correspond to better sensitivity. Top: pixel to pixel sensitivity func- 
tion; regions with bright objects are set to zero. Bottom: smoothed sen- 
sitivity function to show the effects of the dithered single exposures. The 
5<r sensitivity varies between 6.8 x 10~ 21 W m~ 2 pix _1 (white) to 
(black). The contours shown are at 0.9, 1.4, 



3.1 X 10"^ u W m-^pix 
and 1.9X" 20 Wm" 2 pix- 



Lya line at this resolution. As the lines may well appear narrower 
due to the noise, we select peaks with a maximum value over 4a 
above the previously determined fit to the continuum level and a 
2-pixel integrated flux greater than 5<r. Taking a lower limit than 
4a for the peak value of the line leads to too much confusion with 
peaks in the noise; furthermore the detection completeness of the 
low-flux lines would be very low. The integrated 5a flux limit will 
result in the detection of emission lines with a lower limit to the 
flux varying over the whole field between l.OxlO - W m -2 and 
4.4xl0~ 20 W m -2 for an unresolved emission line, as derived 
from the sensitivity function. We only investigate the wavelength 
ranges 4050A < A < 5520A, 5625A < A < 5850A, 5950A 
< A < 6130A, and 6550A < A < 6800A to avoid skylines and 
distortion by bad flux calibration. After selecting all peaks in the 
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Table 1. Observed Lya emitters with their observational properties. The typical errors on the quoted positions are 0.3" and reflect the size of the IFU-fibres. 
The wavelength, flux, and their errors are calculated by a repeated Gaussian fit to the emission lines as described in section |4**T1 A star symbol (*) following 
the candidate number signifies an extended source (i.e. the emission line is spread over more than one IFU fibre). All luminosities Lii nc quoted assume that 
the line lies at the redshift z, given in column 6. 



Nr. 


RA (J2000) 


Dec 


A 




z 






h m s 


o / // 


A 


10" 20 Wm" 2 




10 35 w 


532* 


9:45:30.64 


-24:28:26.9 


5752±1 


3.0±0.7 


3.7315±0.0008 


4.1 ±0.9 


1315* 


9:45:30.95 


-24:28:18.2 


5140±1 


6.7±1.8 


3.2281±0.0008 


6.4±1.7 


1702* 


9:45:31.13 


-24:28:53.7 


5699±1 


2.7±0.8 


3.6879±O.O0O8 


3.6±1.1 


2919 


9:45:31.62 


-24:28:49.0 


5437±1 


1.4±0.4 


3.4724±0.0008 


1.6±0.5 


3210* 


9:45:31.76 


-24:29:15.1 


4676±1 


5.5±1.0 


2.8464±0.0008 


3.9±0.7 


3510 


9:45:31.85 


-24:28:21.5 


6601±1 


2.1±0.5 


4.4299±0.0008 


4.3±1.0 


5644* 


9:45:32.74 


-24:29:05.8 


4766±1 


7.3±1.4 


2.9205±0.0008 


5.5±1.1 


6164 


9:45:32.96 


-24:29:25.9 


4341±1 


3.9±0.8 


2.5709±0.0008 


2.1 ±0.4 


6537 


9:45:33.10 


-24:28:57.1 


6619±1 


1.5±0.3 


4.4447±0.0008 


3.1±0.6 


8060 


9:45:33.72 


-24:29:08.4 


4861±2 


1.7±0.6 


2.9986±0.0016 


1.4±0.5 


8599 


9:45:33.95 


-24:29:15.8 


4280±1 


1.8±0.5 


2.5207±0.0008 


0.9±0.3 


8827 


9:45:34.03 


-24:29:10.5 


4234±2 


2.1±0.7 


2.4829±0.0016 


1.1 ±0.4 


8893 


9:45:34.03 


-24:28:26.2 


461 l±l 


2.3±0.5 


2.7930±0.0008 


1.5±0.3 


10913 


9:45:34.88 


-24:29:13.1 


4232±1 


2.3±0.6 


2.4812±0.0008 


1.2±0.3 



spectra that satisfy the imposed criteria, the peaks that are clearly 
caused by bad pixels, poor flux calibration or other irregularities in 
the data are discarded. We verify the selection procedure by per- 
forming a number of tests. First, the same selection process is run 
again with a smoothed sensitivity function (see bottom Fig. |2j in- 
stead of the pixel-to-pixel values (top Fig.[2j. In this way we check 
if any candidates are selected or rejected due to errors in the noise 
determination in individual fibres; we also run this test without fit- 
ting the continuum, to make sure there are no line candidates se- 
lected or discarded through an erroneous fit. The spectra of the re- 
maining line candidates are then inspected by eye. To this point the 
RMS of the least noisy parts of the spectrum have been used to 
identify candidate line emitters. However, some parts of the spec- 
trum are noisier, therefore the local RMS of the noise around the 
line candidate is determined to check if the peak value is still above 
4a. We then calculate the integrated flux over the two highest pix- 
els to ensure that it is more than 5 times the RMS of the integrated 
local noise. 

Bad pixels and cosmic rays affect the final image even in the 
data cube consisting of the combined frames. Moreover first order 
superpositions are not well removed during the data reduction nor 
are the reduced frames always clear of sky line residues. Therefore 
we construct two more data cubes, comprising 8 and 10 original 
exposures respectively, and calculate at what level the line should 
be detected if it is real. If the line is clearly present in one cube, but 
not in the other, it is likely to be an artefact seen in only one or a few 
original frames. Subsequently for each line candidate all the single 
exposures are visually inspected for artefacts. The final check is to 
see if the line could be the result of crosstalk between fibres, which 
is caused by overflowing of the light from a fibre into a neighbour- 
ing fibre on the CCD. These neighbouring fibres do not necessarily 
correspond to adjoining positions on the sky as there is a compli- 
cated mapping of fibres from the IFU-head onto the CCD. If a line 
candidate passes all these tests, it is declared to be a real emission 
line. After all of these consistency checks we find 15 emission-line 
galaxies. 



4 RESULTS 

4.1 Properties of the line emitters 

Our IFU data facilitates examination of both the spectral and spatial 
properties of the 15 line emitters we have found. The spectra of the 
14 candidate Lya emitters are shown in Fig. [3] The other source is 
a narrow-line active galaxy which is discussed fully in Jarvis et al. 
(2005). 

For each candidate Lya emitter we measure the total flux of 
the emission line by integrating over the spectra of all fibres sam- 
pling the source. The flux of each emission line is measured by 
fitting a Gaussian profile to the emission line. To determine the 
uncertainty in the line profile, we recalculate our fit 1000 times, 
each time adding to the spectrum random noise with an RMS cor- 
responding to the local noise in the spectrum at the position of the 
emission line. The flux and the wavelength of the line are taken 
to be the mean of the 1000 calculations. The standard deviations 
correspond to the la uncertainties on the flux and wavelength. The 
results of these fits can be found in TableQ 

With broad-band images of this field 1 , we derive the contin- 
uum level of the candidate line emitters by measuring (or placing a 
lower limit on) their broad-band magnitude. The sources are taken 
to be detected if the flux within the IFU-fibre aperture is larger 
than 2a above the background level. The V-band image is shown 
in Fig.|4j the locations of the line-emitters are indicated by circles 
which reflect the errors on the positions as quoted in Table Q To 
verify the correctness of the identification of line emitters with ob- 
jects in the broad-band images we compare the coordinates of six 
bright objects that are easily visible in the IFU data cube as well as 
the broad-band image. The mean difference between the two posi- 
tions for each object is ~ 0.3" in RA. In declination however the 



1 These are (1) a B-band image with a 2a limiting magnitude of 
m B limit = 28.4, taken with the Fors2 instrument on the VLT, (2) a V- 
band image with m^/; imit = 28.2, taken with the LRIS camera on the 
Keck Telescope and (3) an I-band image with mj limit = 25.1, also taken 
with LRIS on the Keck Telescope. 
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Figure 3. Spectra of the identified line emitters. The positions of the emission lines are denoted by an arrow. The regions that were excluded from our analysis 
are shaded; at these wavelengths skylines or bad flux calibration dominate the spectra. The redshifts all assume that the emission line is Lyct . 



error varies from ~ 0.3" in the centre of the field to ~ 1" at the 
edges, and we return to this point in section l42l 



4.2 Identification of the emission lines 

It is possible not all our objects are actually Lya emitters at high 
redshift, such that some may be low-redshift interlopers. Within 
our wavelength range any emission line is most likely to be either 
Lya A1216, [Oil] A3727, [OIII] A5007, or Ha A6563. For each 
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Figure 3. cont. 



spectrum we inspect the plausibility of these four cases. The red- 
shift of the object is calculated for each possible line and subse- 
quently the redshifted wavelength of other emission lines are in- 
spected by checking whether the signal at these positions is greater 
than 2a above the local noise, which would indicate a possible 
detection. None of the spectra were found to exhibit any other 
lines under this criterion, which rules out the [OIII] A5007 and 
Ha A6563 possibilities. However confusion with the [Oil] A3727 
line may still occur, since we would not expect to see other bright 
emission lines in our spectra if this is the observed line. The most 
distinguishing characteristic of the [Oil] A3727 line is the fact that 
it is actually a doublet with lines at A = 3726A and A = 3728A; 
unfortunately due to our low resolution we cannot resolve this fea- 
ture. We investigate how likely our lines are to be [Oil] by calculat- 
ing how many [OH] emitters we would expect to find in our volume 
according to Hogg et al. (1998), who carried out a survey at z ~ 1. 
For our observed wavelength range the [Oil] emitters within our 
field of view would be at 0.087 < z < 0.825, corresponding to 
a volume of 214 Mpc 3 (here we use £Im = 0.3, = 0, and 



Ho — 100 km s _1 Mpc~ 1 to facilitate comparison with the [Oil] 
survey of Hogg et al., 1998). Using the [Oil] luminosity function 
combined with the sensitivity function derived from our data, we 
expect to detect 3 [Oil] emitters in our volume. Indeed we have 
found one [Oil] emitter at 2 = 0.62 (see section |3~T1 among the 
bright extended objects that we have excluded from our selection 
procedure. This galaxy could easily be identified as such by the 
clear presence of the 4000A break in its continuum. 

As derived from stellar synthesis models a Lya line, origi- 
nating from a young star-forming regi on, has a typical rest- frame 
equivalent width (EW ) of 100-200A, ICharlot & Failll993h . Sur- 
veys for Lyman Break Galaxies however of ten observe far lower 
equivalent widths for the Lya line. At 2 = 3 lShaplev et alji2003l) 
find that only half of the Lyman Break Galaxies with observed 
Lya emission show an EWo of > 20A and very few objects 
have EWo > IOOA. Also at 2 = 6 Lyma n Break Galaxies usu - 
ally show a detected EW of ~ 20 - 30A fetanwav et alll2004) . 
Nevertheless observations of Lya emitters at high redshift show 
that many of these objects show a Lya EWo of > 200A (e.g. 
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09" 45'" 36.0 s 35.0 s 34.0 s 33.0 s 

Right Ascension 
Figure 4. Keck V-band image of our field of view. The positions of the line 
emitters are marked by a circle the size of an IFU fibre. The accuracy in RA 
is ~ 0.3"; in Dec the error increases with distance from the centre from 
0.3" to up to ~ f". This compromises the certainty of the identification of 
object 3510, 8599, and 8827. 



iMalhotra & Rhoadsll2002l) . Studies of [Oil] emitters on the other 
hand sho w that their equivale nt widths are considerably lower. For 
example Ijansen et alj 120001) have carried out a survey of nearby 
field galaxies, which resulted in equivalent width measurements 
with a maximum of ~ 80 A for [Oil]. Measuring the EW of our 
emission lines thus provides a crucial test of whether the line is 
likely to be Lya or [Oil]. Using the flux of the lines and the (up- 
per limit on) the continuum level, we calculate the (lower limit on) 
equivalent widths of the lines; these are shown in Table [2] For the 
objects detected in the band in which their emission line falls, we 
can use the broad-band images to check whether we can observe a 
spectral break due to hydrogen in the IGM. This is impossible for 
the objects with emission lines in the B-band as we have no bands 
at shorter wavelength. We have however one object (no. 1315) with 
an emission line in V-band t hat is also detecte d in that band; this 
object has (B - V) > 1.3. lDickinsor] ll998t) imposes a limit of 
(B — V) > 1.5 to select Lyman Break Galaxies, which means that 
our B-band image is not deep enough to conclusively determine 
if we see a strong spectral break, although the lower limit to the 
object's colour is close to the selection criterion. The broad-band 
identification for objects 3510, 8599, and 8827 are insecure due to 
their proximity to faint objects in the broad-band images in regions 
of the field where our astrometry becomes less secure ( A5 ~ 1 arc- 
sec). Therefore we have calculated their EWs assuming both no 
identification and an identification with the source found in close 
proximity in the broad-band image. 

As shown in Table [2] all but two of our unknown emission 
lines would have rest-frame equivalent widths of several hundred 
A if they we r e [Oil ], The low-redshift survey for [Oil] -emitters by 
Ijansen et alj l2000l) did not yield any [Oil] emission li nes with an 
EWp l arger than 80A. In the survey at redshift z ~ 1 bv lHogg et alj 
j 19981) the highest detected EW of [Oil] is still below 1 lOA. It is 
therefore highly improbable that any of our high equivalent width 
emission lines are [Oil]. Nevertheless there are two emission lines, 
in objects 3510 and 6537, which could plausibly be [Oil], as the 
lower limit on their EWs would be ~ 119A and ~ 84A respec- 



Table 2. The equivalent widths of the emission lines of our candidate 
Lyo emitters. All values are lower limits unless preceded by '=', which 
represents the actual EW instead of the lower limit. EW D (, S is the observed 
equivalent width, EWq^oii] is rest-frame EW assuming the observed line 
is [Oil] and EWg^Lya is the rest-frame EW if the line is Lya. Quoted errors 
are the result of the uncertainty in the measured flux, magnitude, and red- 
shift. Objects 3510, 8599, and 8827 have an insecure identification in the 
broad-band images. For object 3510 this makes no difference for the calcu- 
lation of the EW as the nearby source is not detected in I-band (in which 
wavelength range its emission line falls). For objects 8599 and 8827 the EW 
is shown for both the case where the nearby source is (1) unassociated and 
(2) associated with the line emission. 



Nr. 


EW ots 


EWo.on 


EWo,Lya 




io 2 A 


io 2 A 


10 2 A 


532 


11±3 


7±2 


2.3±0.5 


1315 


= 13±5 


= 9±4 


= 3.0±1.2 


1702 


10±3 


6±2 


2.0±0.6 


2919 


7±2 


5±1 


1.7±0.5 


3210 


12±2 


10±2 


3.2±0.6 


3510 


2.1±0.5 


1.2±0.3 


0.39±0.09 


5644 


= 10±5 


= 8±4 


= 2.6±1.2 


6164 


16±3 


14±3 


4.6±0.9 


6537 


1.5±0.3 


0.8±0.2 


0.28±0.05 


8060 


7±3 


5±2 


1.8±0.6 


8599 1 


8±2 


7±2 


2.1±0.6 


8599 2 


= 2.9±0.6 


= 2.5±0.5 


= 0.8±0.2 


8827 1 


9±3 


8±3 


2.5±0.8 


8827 2 


= 1.9±0.3 


= 1.7±0.3 


= 0.56±0.09 


8893 


10±2 


8±2 


2.5±0.6 


10913 


10±3 


8±2 


2.8±0.7 



tively which is high but not impossibly so. However it needs to be 
noted that this does not stem from a faint line flux, but from the 
I-band magnitude limit which is much brighter than those of the 
other broad-band images. This results in a lower EW limit. More- 
over the objects show no continuum or any other features in their 
spectra, and it is therefore likely that these objects are Lya emit- 
ters too. Thus for the remainder of this paper we will assume that 
all of our selected candidates are indeed Lya emitters. The inferred 
redshift and line luminosities for our objects are shown in TableQ 



5 DISCUSSION AND CONCLUSIONS 

5.1 The Lya emitter luminosity function 

Our field of view of 1.2' x 1.2' gives a solid angle of 1.5 square 
arcmin. However, for our high-redshift Lya search we have omit- 
ted the regions which contain bright objects. Accounting for this 
reduction in area due to the bright sources, we are left with a solid 
angle of 1.36 square arcmin. If we assume we are only looking at 
Lya emitters, our observed wavelength ranges (see section l3~2l cor- 
respond to redshift ranges of 2.33 < z < 3.54, 3.63 < z < 3.81, 
3.89 < z < 4.04, and 4.39 < 4.59, which yield a total comoving 
volume of 7.8 x 10 3 Mpc 3 . 

The redshift range is large, thus the average flux limit over our 
field of view of 1.8 x 10~ 20 W m -2 (see Fig.|2j corresponds to a 
luminosity limit which rises significantly throughout the depth of 
our volume. If we calculate the average luminosity limit over our 
field we find that it varies within the extremes of our redshift range 
from 8 x 10 34 Watz = 2.3, to 4 x 10 35 Watz = 4.6. 
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Figur e 5. The numb er density of Lya em itters plotted against the luminosity. The filled symbols mark surveys with an average redshift simil ar to ours 
JCowie£^ftJ_ 1998t) : triangles, iKudritzki et alj l2OO0t) : circles) and the open symbols stand for surveys at redshift 2 = 5.7 lAiiki et a?] 120031) : squares, 
iHu eHuTl2004l) : upside-down triangles). Overplotted are two Schechter luminosity functions: the solid line is the fit to all our data points and the dashed 
line is the fit to our two highest luminosity data points and those of the surveys at similar redshift with L > 5 X 10 35 W (dashed horizontal line) to ensure 
completeness. The dotted horizontal lines mark the detection limits of the surveys; for Ajiki et al. (2003) this is the completeness limit. 



Due to the unique nature of our observations we are able to 
probe the Lya luminosity function across a wide range in luminos- 
ity, albeit with a broad redshift bin. In Fig.[3]we show the Lya lumi- 
n osity function derived f rom our observations togeth er with those 
oflKudritzki et alJfeOOOl hereafter KuOO) at z sa 3.1.|Cowie & HvJ 
\ 19981 hereafte 7cH98) at z = 3.4, and lAiiki et alj <2003l her" 
after Aj03) and lHuetalJ <2004 hereafter Hu04) at z » 5.7. We 
have determined the volume of each of the luminosity bins by cal- 
culating the volume within which a source of that luminosity can 
be observed with our flux limit, which varies over the field of view 
according to the sensitivity function derived in section |3~T1 Fig. [5] 
shows that we probe down to much lower luminosities than the 
narrow-band surveys at z > 3; this stems from the fact that at the 
lower redshift limit of our volume we have a significantly lower lu- 
minosity limit. The lack of high-luminosity sources in our volume 
is caused by the small size of the area surveyed. 

To determine if our luminosity function differs from those 
found by the previous surveys at both similar redshift and higher 
redshift, we fit a Schechter function to the data: 

*(^ = **(f) Vv "<f)' 

where is the comoving space density of sources per 

A log Z/Lya, $* is the normalisation, L* is the characteristic break 
luminosity, and a is the slope of the luminosity function. As our 



binned data-points are few, we choose not to fit the Schechter func- 
tion with three free parameters. Therefore we assume a = —1.60: 
the slope of the Schech ter function for z ~ 3 Lyman Break Galax- 
ies ISteidel et allll999h . which also fits well to the luminosity dis- 
tribution of 2 ~ 3 Lya emitters ISteidel et alll200fj) . We determine 
the best fit by carrying out a least-squares minimisation using L* 
and $* as free parameters. The best fit to our data is given by L* = 
5.0 ± 1.8 x 10 35 W and $* = 0.0012 ± 0.0005 Mpc~ 3 (shown 
in Fig. |5|by the solid curve). To ensure completeness we also fit 
a line to a dataset comprising our two highest-luminosity points 
and the results from CH98 and KuOO at L > 5 x 10 35 W, which 
are at similar redshift. We then obtain L* = 5.4 ± 1.7 x 10 35 W 
and $* = 0.0022 ± 0.0012 Mpc" 3 , shown by the dashed curve 
in Fig. [5] The two fits therefore agree with each other within the 
uncertainties. 

It appears that the results of the higher-redshift surveys are 
largely consistent with the lower-redshift Lya luminosity func- 
tion, although the lowest-luminosity results are hampered by in- 
completeness. AJ03 have inferred the same fact by comparing their 
luminosity function with the results of CH98. Hu04 also report that 
the properties of their z = 5.7 sample s hows the same properties as 
their sample at z m 3 l Hu et allll99^) . The inferred conclusion is 
that the luminosity function of L ya emitters does not significantly 
change from z w 3.4 to z = 5.7. lMalhotra & Rhoadsl j2004l) have 
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considered the evolution of the Lya luminosity function at even 
higher redshift, comparing samples at z = 5.7 and z = 6.5. They 
conclude the luminosity functions do not statistically deviate from 
each other, implying that the reionisation of the Universe must have 
been largely complete at z = 6.5. 

We can calculate how many high-luminosity sources we 
would expect to find according to our luminosity function. In this 
way we can check our earlier statement that we find none of these 
objects because our probed volume is too small. If we integrate the 
luminosity function from the edge of our highest luminosity bin 
to infinity over the volume probed by our survey, we find that we 
would expect to find only one object with a higher Lya luminos- 
ity. The lack of high-luminosity sources in our survey is therefore 
within the Poisson error of our expectation and is not significant. 



5.2 The star-formation rate 

As the Lya line is excited when interstellar gas is photoionised by 
the UV radiation of young and hot blue stars, the Lya luminos- 
ity of a galaxy is directly related to its star-formation rate (SFR). 
We can derive this relatio n by using the case B recombination the- 
ory iBrocklehursl Il97ll) to deduce the Ha luminosity of a star- 
fo rming galaxy from its Lya luminosity, and the equation found 
bv lKennicuttlll998h to link in turn its star-formation rate to its Ha 
luminosity. The formulae are respectively Z/Lya = 8.7Lu a , and 
SFR = 7.9 x 10" 35 L Ha M yr~\ where L Hq is the Ha lumi- 
nosity in Watts and continuous star formation is assumed with a 
Salpeter Initial Mass Function from 0.1 — 100 Mq. Thus we get: 

SFR = 9.1 x 10" 36 f^ 1 ) M yr-\ 

Using this relation we estimate the star-formation rate for each 
of our objects. It is important to note that these star-formation rates 
are calculated under the assumption that the Lya line is solely 
photoionised by hot stars and neglecting any extinction. The SFRs 
range from about 1 to 6 solar masses per year, with an average of 
2.6 ±1.6 M yr _1 . 

Summing the individual SFRs of our objects and dividing by 
their corresponding volume gives a comoving star-formation rate 
density of psfr - 0.0067 ± 0.0005 M yr _1 Mpc -3 . Integrating 
over the luminosity function yields Psfr = 0.008 Moyr^Mpc- 3 
for L* = 5.0 x 10 35 W and $* = 0.0012 Mpc~ 3 , which ap- 
plies to our own dataset, and psfr ~ 0.017 M0yr~ 1 Mpc~ 3 for 
L* = 5.4 x 10 35 W and $* = 0.0022 Mpc~ 3 as derived from the 
combination of our dataset together with CH98 and KuOO. We com- 
pare this SFR densit y to those from previous studies, compiled by 
iKodaira et alj 1200 3l) . recalculated with our cosmological parame- 
ters and more recent surveys (Fig. [6] Table[3j. 

Although our result is lower than those of multi-colour sur- 
veys, it is in good agreement with the Lya surveys carried out 
around the same redshift. The discrepancy between the Lya 
searches and the multi-colour surveys is probably because the Lya 
searches only set a lower limit to the SFR density. This is due to the 
fact that the SFR based on Lya luminosity does not take account 
of extinction by dust and interstellar and intergalactic gas. This is 
a very important effect as Lya is a resonant line and the luminos- 
ity of the Lya line can be severely affected by absorption from 
HI. However, correcting for obscuration requires extensive knowl- 
edge of the amount of dust and gas in the light path and is there- 
fore very uncertain and deep optical and near-infrared observations 
are needed to make these corrections. To compare our result with 



the multi-colour surveys, we can apply a correction for the differ- 
ence between the UV continuum emission and the Lya emission, 
which alleviates some of the problems associated with HI absorp- 
tion of the Lya line . We use SFR(Lya) = 0.77 x SFR(UV) 
(see lHuetalll2002l) . where SFR(UV) is the star-formation rate 
estimate derived from the UV continuum emission. The corrected 
SFR density then becomes psfr = 0.010 Meyr -1 Mpc~ 3 for our 
data and psfr = 0.022 Moyr -1 Mpc~ 3 for the combined dataset. 
These agree well with the SFR densities found by multi-colour sur- 
veys at similar redshift, although we note the dust obscuration may 
still be a severe problem. 

The lack of any evidence for evolution of the luminosity func- 
tion of Lya emitters between z w 3.4 and z = 5.7 suggests that the 
star-formation rate density does not change in this redshift interval, 
although our survey is severely limited by low number statistics. In- 
spection of Fig.|6|reveals that various surveys around z — 6 show 
evidence for a decrease in the cosmic star-formation rate density to- 
ward high redshift. For Lya-emitter surveys this could be caused by 
the increasing optical depth of HI to the Lya emiss i on line. How- 
ever, the multi-colour surveys IStanwav et allEoO^lBunker et all 
|2004) which calculate the cosmic star-formation rate density by 
integrating over the luminosity function of Lyman Break Galax- 
ies, also point to a decrease in the star-formation density. A future 
wider-field IFU survey would result in a large sample of Lya emit- 
ters discovered in a perfectly consistent manner. This would facil- 
itate investigation of the star-formation rate density within a large 
redshift range by comparing exactly the same sort of objects with 
the same selection criteria. This could provide significant clues to 
the evolution of the star-formation rate density throughout the Uni- 
verse. 

5.3 Clustering around the radio galaxy 

It is interesting to determine whether the fact that our observations 
were centred on radio galaxy MRC 0943-242 at z = 2.92 has any 
influen ce on the number of Lya emitters we find. IVenemans et alj 
120021) suggest that luminous radio galaxies are tracers of regions 
of galaxy overdensities in the early Universe. They find Lya emit- 
ters around the radio galaxy TN 11338-1942 at z = 4.1 with 
a relative overdensity of ~ 15 compared to blank-field surveys. 
iKudritzki etafl <2000l) carried out a Lya search with a similar flux 
limit to our average limit, targeted at a blank field at z — 3.1 which 
is close to the redshift of our radio galaxy. They derive a number 
density of 14400Az~ 1 deg~ . Based on this we would expect to 
find ~ 1 Lya emitter at redshift z ~ 2.92 over our entire field, 
within a redshift interval of Az ~ 0.004. Indeed we find one emit- 
ter with 2 rg - z em = 0.002 (object 5644). 

It is evident that the angular size of our field is much too small 
to determine whether there is an overdensity around our target. Al- 
though it shows that we can compare our results with blank field 
studies without a significant bias due to the presence of the radio 
galaxy. 

We also note that three of our Lya emitters lie withing Az = 
0.04 of each other at z ~ 2.5. It is difficulty to say that this is def- 
initely an overdensity, but it is suggestive and further observations 
are needed. However, this shows the power that IFU-surveys may 
have in identifying regions of galaxy overdensities at high redshift. 



6 SUMMARY 

We summarise our results and conclusions as follows: 
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Figure 6. Star-formation rate densities as derived by various types of surveys. The result from this work derived from integrating over the luminosity function 
fit to our data alone is denoted by the star symbol. The different types of surveys are marked with different symbols: the open circles are Ha searches, the 
open triangles are surveys aimed at oxygen emission lines, the filled triangles are multi-colour surveys, and the open squares are Lya searches. All the SFR 
densities and their relevant references can be found in Tablel3l The x-error bars reflect the redshift range of the surveys. 



Table 3. The star-formation rate densities found by previous studies, shown in Fig.lol 



Label 


z 


PSFR 




Survey Type 


Source 


Label 


z 


PSFR 




Survey Type 


Source 






M Q yr~ 


l 










M e yr~ 


i 






A 


0.02 


1.2x10 


-2 


Ha 


Galleso et al. (19961 


P 


2.75 


2.9x10" 


-2 


Multi-colour 


Madau et al. (19981 


B 


0.15 


2.5x10 


-2 


Ha 


Tresse & Maddox ( 19981 


Q 


3.04 


5.0x10" 


-2 


Multi-colour 


Steidel et al. (19991 


C 


0.25 


2.4x10 


-2 


Ha 


Hippelein et al. (20031 


R 


3.14 


6.4x10" 


-3 


Lya 


Kudritzki et al. (20001 


D 


0.25 


1.2x10 


-2 


Multi-colour 


Treveret al. (19981 


S 


3.44 


4.7x10" 


-3 


Lya 


Cowie & Hu (19981 


E 


0.35 


1.1x10 


-2 


Multi-colour 


Lilly etal. (19961 


T 


3.45 


8.0x10" 


-3 


Lya 


This paper 


F 


0.40 


2.4x10 


-2 


[OIII] 


Hippelein et al. (20031 


U 


3.72 


4.1x10" 


-4 


Lya 


Fuiita et al. (20031 


G 


0.63 


2.2x10 


-2 


Multi-colour 


Lilly et al. (19961 


V 


4.00 


1.1x10" 


-2 


Multi-colour 


Madau et al. (19981 


H 


0.64 


7.2x10 


-2 


[OIII] 


Hippelein et al. (20031 


w 


4.13 


4.0x10" 


-2 


Multi-colour 


Steidel et al. (19991 


I 


0.75 


4.2x10 


-2 


Multi-colour 


Connolly et al. (19971 


X 


4.85 


3.5x10" 


-2 


Multi-colour 


Iwata et al. (20031 


J 


0.88 


4.2x10 


-2 


Multi-colour 


Lillvetal. (19961 


Y 


4.86 


6.3x10" 


-3 


Lya 


Ouchi et al. (20031 


K 


0.88 


1.1x10 


-1 


[Oil] 


Hippelein et al. (20031 


Z 


5.70 


1.2x10" 


-3 


Lya 


Aiiki et al. (20031 


L 


1.20 


2.3x10 


-1 


[on] 


Hippelein et al. (20031 


AA 


5.74 


5.0x10" 


-4 


Lya 


Rhoads et al. (20031 


M 


1.25 


5.8x10 


-2 


Multi-colour 


Connolly et al. (19971 


BB 


5.80 


0.3x10" 


-2 


Multi-colour 


Stanwav et al. (20031 


N 


1.75 


4.4x10 


-2 


Multi-colour 


Connolly etal. (19971 


CC 


6.00 


0.5x10" 


-2 


Multi-colour 


Bunker et al. (20041 


O 


2.21 


1.2x10" 


-1 


Ha 


Moorwood et al. (20001 


DD 


6.57 


5.2x10" 


-4 


Lya 


Kodaira et al. (20031 



• We have discovered 14 Lya emitters with a flux varying from 
1.4 to 7.3 x 1(T 20 W irT 2 in a comoving volume of 7.8 x 10 3 Mpc 3 . 
The resulting number density is 0.0018±g;oSoi Mpc~ 3 . 

• Fitting a Schechter luminosity function to the data shows that 
the luminosity distribution of the sources is comparable to that of 



Lya emitters at z — 5.7. This implies that we do not see any evo- 
lution of Lya emitters between z ~ 3.4 and z = 5.7. 

• The star-formation rate of the individual objects lies within 
a range from about 1 to 6 solar masses per year. The derived 
cosmic star formation rate density is psfr = 6.7 ± 0.5 x 
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10~ 3 M0yr -1 Mpc -3 . This is consistent with the results of Lya 
emitter searches at similar redshift. It is however significantly lower 
than the results of multi-colour surveys. This is due to the fact 
that the Lya searches only yield a lower limit to the cosmic star- 
formation rate density because of absorption of the Lya emis- 
sion and the lack of a reliable luminosity function. We can com- 
pare our result to the multi-colour surveys if we determine the 
cosmic star-formation rate by integrating our fit of the Schechter 
luminosity function to the combination of our dataset with sur- 
veys at similar redshift and correcting for the difference of extinc- 
tion between Lya luminosity and UV luminosity. We then obtain 
Psfr ~ 2.2 x 10~ 2 M Q yr~ 1 Mpc -3 , which is similar to multi- 
colour results. 

• We find one companion to the central radio galaxy at z = 
2.92. We do not find further evidence for clustering around the ra- 
dio galaxy, as our field of view is too small to allow for such investi- 
gations. However, we do find three Lya emitters within Az = 0.04 
at z ~ 2.5 which suggests that there may be an overdensity of 
galaxies at this redshift, although further observations are needed 
to confirm this. 

• An important conclusion to this project is that the IFU- 
technique is well suited to detect Lya emitters and trace the star- 
formation history of the Universe. It is however vital to carry out 
larger surveys to probe a larger volume and therefore detect a higher 
number of sources. It will then become possible to reliably sample 
the luminosity density of Lya emitters throughout a large redshift 
range with a single survey. This will facilitate studying the evolu- 
tion of the properties of Lya emitters throughout time and therefore 
map the evolution of the cosmic star-formation rate density reli- 
ably, through comparison of galaxies selected by exactly the same 
method. 



ACKNOWLEDGMENTS 

We thank Wil van Breugel for the B-, V- and I-band images and the 
LRIS spectrum. We would also like to thank Steve Rawlings, Huub 
Rbttgering and Richard Wilman for useful discussions. The IFU 
data published in this paper have been reduced using VIPGI, devel- 
oped by INAF Milano, in the framework of the VIRMOS Consor- 
tium activities. We thank Bianca Garilli and Marco Scodeggio for 
their help regarding VIPGI. CVB would like to acknowledge the 
financial support from the ERASMUS program, the International 
Study Fund of Leiden University (LISF) and a PPARC studentship, 
MJJ acknowledges the support of a PPARC PDRA. 



References 

Ajiki M., Taniguchi Y., Fujita S. S., Shioya Y., Nagao T, Mu- 
rayama T, Yamada S., Umeda K., et al., 2003, AJ, 126, 2091 

Becker R. H., Fan X., White R. L., Strauss M. A., Narayanan 
V. K., Lupton R. H., Gunn J. E., Annis J., et al, 2001, AJ, 122, 
2850 

Brocklehurst M., 1971, MNRAS, 153, 471 

Bunker A. J., Stanway E. R., Ellis R. S., McMahon R. G, 2004, 

MNRAS, 355, 374 
Chariot S., Fall S. M., 1993, ApJ, 415, 580 
Connolly A. J., Szalay A. S., Dickinson M., Subbarao M. U., 

Brunner R. J., 1997, ApJ, 486, LI 1 
Cowie L. L., Hu E. M., 1998, AJ, 115, 1319 



Dickinson M., 1998, in The Hubble Deep Field Color-Selected 

High Redshift Galaxies and the HDF. p. 219 
Djorgovski S. G, Castro S., Stern D., Mahabal A. A., 2001, ApJ, 

560, L5 

Eggen O. J., Lynden-Bell D., Sandage A. R., 1962, ApJ, 136, 748 
Fujita S. S., Ajiki M., Shioya Y, Nagao T, Murayama T, 

Taniguchi Y, Okamura S., Ouchi M., et al., 2003, A J, 125, 13 
Gallego J., Zamorano J., Aragon-Salamanca A., Rego M., 1996, 

ApJ, 459, L43 
Gunn J. E., Peterson B. A., 1965, ApJ, 142, 1633 
Hippelein H, Maier C, Meisenheimer K., Wolf C, Fried J. W., 

von Kuhlmann B., Kiimmel M., Phleps S., et al., 2003, A&A, 

402, 65 

Hogg D. W., Cohen J. G, Blandford R., Pahre M. A., 1998, ApJ, 
504, 622 

Hu E. M., McMahon R. G, Egami E., 1997, in The Hubble Space 
Telescope and the High Redshift Universe Galaxies at z > 4.5. 
p. 91 

Hu E. M., Cowie L. L., McMahon R. G, 1998, ApJ, 502, L99 
Hu E. M., Cowie L. L., McMahon R. G, Capak P., Iwamuro F, 

Kneib J.-R, Maihara T, Motohara K., 2002, ApJ, 568, L75 
Hu E. M., Cowie L. L., Capak P., McMahon R. G, Hayashino T, 

Komiyama Y, 2004, AJ, 127, 563 
Iwata I., Ohta K., Tamura N., Ando M., Wada S., Watanabe C, 

Akiyama M., Aoki K., 2003, PASJ, 55, 415 
Jansen R. A., Fabricant D., Franx M., Caldwell N., 2000, ApJS, 

126, 331 

Jarvis M.J., van Breukelen C, & Wilman R.J., 2005, preprint 

lastro-ph/0412600! 
Kennicutt R. C, 1998, ARA&A, 36, 189 

Kodaira K., Taniguchi Y, Kashikawa N., Kaifu N., Ando H., 
Karoji H., Ajiki M., Akiyama M., et al., 2003, PASJ, 55, L17 

Koo D. C, Kron R. T, 1980, PASP, 92, 537 

Kudritzki R.-P, Mendez R. H, Feldmeier J. J., Ciardullo R., Ja- 
coby G. H, Freeman K. C, Arnaboldi M., Capaccioli M., et al., 
2000, ApJ, 536, 19 

LeFevre O., Saisse M., Mancini D., Brau-Nogue S., Caputi O., 
Castinel L., D'Odorico S., Garilli B., et al., 2003, in lye M„ 
Moorwood A. F. M., eds., Proceedings of the SPIE, Volume 
4841, Instrument Design and Performance for Optical/Infrared 
Ground-based Telescopes, p. 1670 

Lilly S. J., Le Fevre O., Hammer F, Crampton D., 1996, ApJ, 460, 
LI 

Lowenthal J. D., Hogan C. J., Leach R. W., Schmidt G. D., Foltz 

C. B., 1990, ApJ, 357, 3 
Madau P., Pozzetti L., Dickinson M., 1998, ApJ, 498, 106 
Madau P., Haardt F, Rees M. J., 1999, ApJ, 514, 648 
Malhotra S„ Rhoads J. E., 2002, ApJ, 565, L71 
Malhotra S., Rhoads J. E., 2004, ApJ, 617, L5 
Moorwood A. F. M., van der Werf P. P., Cuby J. G, Oliva E., 2000, 

A&A, 362, 9 

Ouchi M., Shimasaku K., Furusawa H, Miyazaki M., Doi M., 
Hamabe M., Hayashino T, Kimura M., et al., 2003, ApJ, 582, 
60 

Peebles P. J. E., 1971, A&A, 11, 377 

Press W. H., Schechter P., 1974, ApJ, 187, 425 

Rhoads J. E., Dey A., Malhotra S., Stem D., Spinrad H., Jannuzi 

B. T, Dawson S., Brown M. J. I., et al., 2003, AJ, 125, 1006 
Sandage A., Freeman K. C, Stokes N. R., 1970, ApJ, 160, 831 
Scodeggio M., Franzetti P., Garilli B., Zanichelli A., Paltani S., 

Maccagni D., 2004, preprint I astro-ph/0409248 I 



12 Caroline van Breukelen, Matt J.Jarvis & Bram P. Venemans 



Shapley A. E., Steidel C. C, Pettini M., Adelberger K. L., 2003, 
ApJ, 588, 65 

Stanway E. R., Bunker A. J., McMahon R. G., 2003, MNRAS, 
342, 439 

Stanway E. R., Bunker A. J., McMahon R. G., Ellis R. S., Treu T., 
McCarthy P. J., 2004, ApJ, 607, 704 

Steidel C. C., Adelberger K. L., Giavalisco M., Dickinson M., Pet- 
tini M., 1999, ApJ, 519, 1 

Steidel C. C, Adelberger K. L., Shapley A. E., Pettini M., Dick- 
inson M., Giavalisco M., 2000, ApJ, 532, 170 

Thompson D. J., Smith J. D., Djorgovski S., Trauger J., 1992, 
Bulletin of the American Astronomical Society, 24, 806 

Tresse L., Maddox S. J., 1998, ApJ, 495, 691 

Treyer M. A., Ellis R. S., Milliard B„ Donas J., Bridges T. J., 
1998, MNRAS, 300, 303 

Venemans B. P., Kurk J. D., Miley G. K., Rottgering H. J. A., van 
Breugel W., Carilli C. L., De Breuck C, Ford H., et al., 2002, 
ApJ, 569, Lll 

White S. D. M., Rees M. J., 1978, MNRAS, 183, 341 



